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Jet-cloud collisions in the jet of the Seyfert galaxy NGC 3079 



o 
o 

(N 
X> 

(X, 

OS 



E. Middelberg 1 *, I. Agudo, 2,3 , A. L. Roy 2 , and T. P. Krichbaum 2 

1 Australia Telescope National Facility, PO Box 76, Epping NSW 1710, Australia 
2 Max-Planck-Insitut fur Radioastronomie, Aufdem Hugel 69, 53121 Bonn, Germany 
^Institute de Astroflsica de Andalua'a, CSIC, Apartado 3004, 18080 Granada, Spain 



oo 

(N 
O 

o 

Oh 
6 



ABSTRACT 

We report the results from a six-year, multi-epoch very long baseline interferomertry mon- 
itoring of the Seyfert galaxy NGC 3079. We have observed NGC 3079 during eight epochs 
between 1999 and 2005 predominantly at 5 GHz, but covering the frequency range of 1 .7 GHz 
to 22 GHz. Using our data and observations going back to 1985, we find that the separation of 
two of the three visible nuclear radio components underwent two decelerations. At the time of 
these decelerations, the flux density of one of the components increased by factors of five and 
two, respectively. We interpret these events as a radio jet component undergoing compression, 
possibly as a result of a collision with ISM material. This interpretation strongly supports the 
existence of jets surrounded by a clumpy medium of dense clouds within the first few parsecs 
from the central engine in NGC 3079. Moreover, based on recently published simulations of 
jet interactions with clumpy media, this scenario is able to explain the nature of two addi- 
tional regions of ageing synchrotron material detected at the lower frequencies as by-products 
of such interactions, and also the origin of the kpc-scale super bubble observed in NGC 3079 
as the result of the spread of the momentum of the jets impeded from propagating freely. The 
generalization of this scenario provides an explanation why jets in Seyfert galaxies are not 
able to propagate to scales of kpc as do jets in radio-loud AGN. 
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1 INTRODUCTION 

Active galactic nuclei (AGN) are conventionally divided into radio- 
loud and radio-quiet objects, depending on their ratio of radio-to- 
optical luminosity. Although radio-quiet objects such as Seyfert 
galaxies have luminosities of up to 1000 times less than radio- 
loud objects, they mos tly exhibit high-brightne ss temperature (non- 
thermal) components ( Middelberg et al. 2004), indicating that the 
radio emission is synchrotron emission. However, it is not clear 
why the generation of radio emission in Seyfert galaxies is dif- 
ferent from that in radio-loud objects such as quasars and ra- 
dio galaxies. In particular, Seyferts mostly do not exhibit long 
and confined jets as are typically observed in radio-loud objects. 
Only in few radio-quiet objects have jet-li ke features been ob - 
served, pro minent examples are NGC 4258 iMivoshi et aill995). 



NGC 4151 dMundell et alj 120031), NGC 3079 dlrwin & Se aquist 



Kondra tko et all 120051) . and NGC 1068 1 19851 : lHaschick et all Il990l : iTrotter et alj 1 19981 : iKondratko et j 



Il988l: iTrotter et alj|l998 
dRovetal.lll998h . 

NGC 3079 is a nearby (15.0 Mpc, Ide Vaucouleurs et al.1 

Il99ll Hn = TS kms^Mpc" 1 1 mas yr" 1 = 0.24 c) LINER 



carried out by ITrotter et al. I dl998l), [S awada-S atoh et alj d200Gt) . 
IKondratko etafl d2005l) and iMiddelberg et al.1 (2005) revealed that 
the separation of A and B was increasing at a speed between 0.11c 
and 0.16 c between 1985 and 2000. On the axis connecting A and 
B, indications of another c omponent called C were already found 
by|l rwin & Seaquistl (fl988), With very sensitive observations it has 
recently been possible to detect a structure that c overs about half 
the p rojected distance of 2pc between A and B (Kondratko et al. 
2005), indicating a jet-like feature. NGC 3079 therefore appears to 
be a promising candidate source to investigate the nature of radio 
emission in Seyferts. 



NGC 3079 also contains one of the brightest known extra- 
galctic water masers (i^est = 22.23508 GHz), w hich has been sub- 
ject to investigation for more than a decade (e.g., Haschick & Baan 

ID 



2005). The distribution of the masers and their velocity dispersion 



(Heckman 1980) or Seyfert 2 dSosa-Brito et alj|200ll) galaxy in 
an edge-on, dusty spiral . In the radio, pilot VLBI observations by 
llrwin & Seaquistl J 19881) resolved the core into two strong com- 
ponents, A and B, separated by 20.2 mas. VLBI observations 
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indicate a thick molecular disk rotating around a massive object, 
presumably a black hole. 

We have observed the radio emission in NGC 3079 with the 
VLB A at 5 GHz, 15 GHz, and 22 GHz, using phase referencing in 
most epochs, to monitor precisely its structural evolution. Single 
epoch images at 1.7 GHz and 2.3 GHz were obtained to probe a 
wider spectral range. 
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2 OBSERVATIONS 

The data presented here were obtained with the VLBA over eight 
epochs between 1999 and 2005, partly supplemented by other tele- 
scopes. The sensitivity, (u, v) coverage, and resolution varied due 
to the wide range in frequency. Also, phase referencing was used 
in the later observations, but not in the 1999/2000 observations. A 
summary of the observations is presented in TableQ] 

At 1.7 GHz and 2.3 GHz, NGC3079 was observed during a 
single epoch only, on 22 September 2002, using phase referencing 
to the nearby calibrator 4C+55.17. The Effelsberg telescope was 
included in the 1.7 GHz observations only, to increase the resolu- 
tion at this low frequency. 

We have observed NGC 3079 at 5 GHz during eight epochs 
between November 1999 and August 2005 with the VLBA, in 
part supplemented by the Effelsberg telescope and a single VLA 
antenna. The observations carried out in 2002-2005 were phase- 
referenced to the nearby source 4C+55.17. We have also obtained 
a 5 GHz VLBI image from a 16-station global VLBI snapshot ob- 
servation performed in 1994 (project GK01 1), but were only able to 
reliably measure the separation of A and B from this observation. 

During three epochs in 1999-2000, we have observed 
NGC 3079 at 15 GHz with the VLBA, the Effelsberg telescope and 
a single VLA antenna. The observations yielded only few detec- 
tions on some baselines shorter than 100 MA during the November 
1999 and April 2000 epochs, but in November 2000, the source was 
detected on most baselines shorter than 120 MA. 

During the three epochs in 2004-2005, we observed the wa- 
ter maser emission in NGC 3079 at j/ rest = 22.23508 GHz and 
the 22 GHz continuum emission. The frequency resolution was 
31.25 kHz, corresponding to a velocity resolution of 0.42 kms -1 . 
The continuum emission was phase-referenced to the brightest 
maser feature. 

2.1 Calibration 

The data were calibrated using standard procedures implemented 
in the Astronomical Image Processing System, AIPS. The visibil- 
ity amplitudes were calibrated using noise-adding radiometry and 
established antenna gains, and amplitude errors of the order of 
5 % arising from sampler threshold variations were corrected us- 
ing the autocorrelation data. Phase rotation due to parallactic an- 
gle changes was corrected, and instrumental delays were calibrated 
using fringe-fitting on a short observation of a strong source. For 
phase calibration of the observations carried out in 1999-2000, we 
fringe-fitted the data using a point source model, with solution 
intervals of around 4min, yielding detections on baselines up to 
100 MA at both 5 GHz and 15 GHz. At 22 GHz, the phase cali- 
bration was carried out using six frequency channels with a total 
bandwidth of 187.5 kHz for self-calibration on the strongest water 
maser feature, at a velocity of 956 km s~ 1 . The maser was detected 
on most baselines, and the phase solutions were then applied to all 
frequency channels. 

2.2 Eearth orientation parameter errors and phase 
calibrator position changes 

The observations in 2004 and 2005 were correlated using predicted 
rather than measured earth orientation parameters, a situation aris- 
ing from a software upgrade at the VLBA correlator. Furthermore, 
the position of the phase calibrator 4C +55.17 used for phase refer- 
encing of the 5 GHz observations in the two 2005 epochs differed 



from the position used in earlier epochs. This happened because 
calibrator positions in the program used to schedule VLBA obser- 
vations, SCHED, were taken from a data base at the Array Opera- 
tions Center at the time of observing, rather than from the schedule, 
and hence were subject to unforseeable changes. However, the ef- 
fects of both problems are predictable and precisely known, and 
were easily rectified. 

We retrieved the final earth orienta- 
tion parameters from the USNO solution at 
|http : // gemini . gsf c . nasa . gov/ solvejsave/usno_f inals . erp| 
and applied the corrections to the data using the AIPS task VL- 
BAEOP. This process has been well tested and the residual 
phase errors after c orrection are of the order of < 5° at 5 GHz 
(Walk er et al. ||2003). The position of the calibrator in the 2002- 
2004 observations, 4C+55.17, was modified using the AIPS task 
CLCOR. The right ascension was decreased by 3.375 mas and 
the declination was decreased by 0.792 mas, to agree with the 
presumably more accurate position used in the 2005 observations, 
RA=09:57:38. 184971, Dec=55:22:57.76924 (J2000). 

2.3 Image analysis 

The relative positions and flux densities of the 5 GHz and 15 GHz 
observations in 1994 and 1999-2000 were measured from self- 
calibrated images (because no phase referencing was carried out), 
tapered to a resolution of 3 x 3 mas 2 as a compromise between 
(u, v) coverage and sensitivity. Absolute positions were not mea- 
sured in these observations. 

The absolute positions of the 5 GHz and 22 GHz observations 
in 2002-2005 were measured from purely phase-referenced images, 
with no self-calibration steps carried out, and tapered to 3 x 3 mas 2 . 
The flux densities at 5 GHz were measured from phase-referenced, 
self-calibrated, and tapered images, to remove residual phase er- 
rors and improve the signal-to-noise ratio. Flux densities at 22 GHz 
were measured from purely phase-referenced images because self- 
calibration turned out to be marginal. In self-calibration, the result- 
ing flux densities of components A and B were found to vary by as 
much as a factor of three when the calibration strategy was changed 
only very little and our models did not reliably converge. We there- 
fore consider the 22 GHz flux densities to be very unreliable and do 
not use them in our analysis. 

We have measured the positions of components A, B, and E 
and their peak flux densities by fitting a paraboloid to a square of 
three by three pixels, using the location of the brightest pixel as 
a starting point. Exploring five images from a reasonable imaging 
parameter space revealed a position uncertainty of 0.15 mas (lcr) 
for the 5 GHz and 15 GHz images, and of 0. 1 mas for the 22 GHz 
images. Hence separations, which include two position measure- 
ments, have errors of \ // 2 times the position errors. The errors of 
our 1994 observations are larger because of the very sparse (u, v) 
coverage. 

Integrated flux densities were measured by integrating over 
the source region. We have estimated the flux-density errors as fol- 
lows: peak flux densities have a 5% (10%) scaling error at fre- 
quencies at and below 5 GHz (15 GHz) plus the image rms added 
in quadrature. Integrated flux densities have the same scaling errors, 
plus the image rms in quadrature, plus an error arising from the in- 
tegration over the source region, which was derived experimentally, 
added in quadrature. For components A and F, the integration er- 
ror alone amounts to 1 mjy because they blend with C or are poorly 
defined, and to 0.3 mjy for components B and E. 

For illustration, we show in Figures [TT5l l.7 GHz and 2.3 GHz 
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suppl. Antennas 


At 


Au 


Phase-ref 




GHz 




min 


MHz 




1994-02-28 


5.0 


Eb/Y /On/Me/No 


~10 


56 


n 


1999-11-20 


1.7,5.0, 15 


Yl/Eb,Yl/Eb,Yl 


108, 144, 144 


32 


n 


2000-03-06 


1.7,5.0, 15 


Yl/Eb,Yl/Eb,Yl 


108, 144, 136 


32 


n 


2000-11-30 


1.7,5.0, 15 


Yl/Eb,Yl/Eb,Yl 


108, 144, 144 


32 


n 


2002-09-22 


1.7,2.3,5.0 


Eb 


114, 107, 112 


64 


y 


2004-06-14 


5.0, 22 




199, 94 


64 


y 


2004-11-03 


5.0, 22 




194, 82 


64 


y 


2005-04-04 


5.0, 22 




195, 107 


64 


y 


2005-08-18 


5.0, 22 




192, 107 


64 


y 



Table 1. Summary of the observations. We give the observing dates, the frequencies observed, antennas supplementing the VLBA (Eb - Effelsberg, Yl - one 
VLA antenna, Y - phased VLA, Me - Medicina, On - Onsala 85, No - Noto), the integration times, the bandwidths, and whether phase referencing was used 
or not. 2-bit sampling was used except 1994, when 1-bit sampling was used. 
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Figure 1. 1.7 GHz, phase-referenced and self-calibrated image of 
NGC 3079, observed on 22 September 2002. Axes are relative to RA 
10h01m57.7906s Dec 55°40'47.1788" (J2000), and self-calibration pro- 
duced no measureable shift in the coordinates from phase referencing. Con- 
tours are drawn at 0.18 mjy beam -1 X 2 N (N = 0, 1, 2, ...), and the 
beam size is 4.26 X 2.99 mas 2 . Only components E and F have been de- 
tected, and a wider region is shown to demonstrate that A, B, and C are not 
visible. 

images from the 22 September 2002 epoch, 5 GHz and 22 GHz im- 
ages from the 18 August 2005 epoch, and a 15 GHz image from 
the 30 November 2000 epoch. Details can be found in the image 
captions. 

Caption to Table 2 Component and image parameters. Note 
that the 22 GHz flux densities are not reliable because self- 
calibration was marginal, see the discussion in Section [2~3l 
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Figure 2. 2.3 GHz, phase-referenced and self-calibrated image of 
NGC 3079, observed on 22 September 2002. Axes are relative to RA 
10h01m57.7906s Dec 55°40'47.1788" (J2000), and self-calibration pro- 
duced no measureable shift in the coordinates from phase referencing. Con- 
tours are at 0.39 mjy beam -1 X 2^ (N = 0, 1, 2, ...), and the beam size 
is 5.75 X 4.10 mas 2 . Components E and F are weaker than at 1.7 GHz, 
and A is weak, but detected. 

3 RESULTS 

3.1 The maser emission 

The purpose of observing the maser emission in this project was to 
phase reference the radio continuum emission to a quasi-stationary 
point (cf. the next section). In our data, the 1 a thermal noise per 
channel is 1 2 m Jy beam ~ 1 , only allowing us to detect the brightest 
features. We therefore do not attempt an astrophysical modelling of 
the maser emission. A spectrum showing all detected lines is shown 
in Figure[6] 



3.2 Position of B with respect to the maser emission 

One of the aims of the 22 GHz observations was to measure the 
positions of A and B with respect to the water maser emission. 
We found that the separation of the brightest maser feature at 
956 km s" 1 from the continuum component B was 6.7 mas ± 
0.14 mas, with a position angle of 82 ° . This is in agreement with 



U. 14 mas, witn a position angle ot 82 . Inis is in agreement witn 
meas urements by iTrotter et alj 1 1998) (6.7mas ± 0.1 mas in PA 



81°), [Sawada-Satoh et alj d2002h (« 6.7 mas in PA 81°), and 



iKondratko etal] J2005h (6.74 mas ± 0.13 mas in PA 80.75° ± 
1.04°), which do not reveal any conclusive motion between Jan- 
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Observations at 2.3 GHz 






















2002-09-22 


1.2±0.1 


19.8±1.0 


2.1±0.2 
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21.0±1.2 


7.2±0.7 


0.13 


6.0 
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Observations at 5 GHz 






















1999-11-20 


24.8±1.2 


12.8±0.7 3.1±0.2 




27.1±1.7 


15.9±1.0 4.1±0.6 




0.097 


3.0 


3.0 


0.0 


2000-03-06 


23.3±1.2 


13.1±0.7 3.1±0.2 




26 0±1 6 


18.7±1.1 3.2±0.6 




0.095 


3.0 


3.0 


0.0 


2000-11-30 


22.9±1.1 


13.2±0.7 3.8±0.2 




23.8±1.6 


18.8±1.1 6.4±0.6 




0.092 


3.0 


3.0 


0.0 


2002-09-22 


21.1±1.1 


14.0±0.7 5.4±0.3 




26.0±1.6 


22.4±1.3 7.9±0.7 


6.9±1.1 


0.089 


3.0 


3.0 


0.0 


2004-06-14 


41.7±2.1 


16.6±0.8 5.4±0.3 




47.0±2.6 


21.8±1.2 8.9±0.7 




0.085 


3.0 


3.0 


0.0 


2004-11-03 


41.7±2.1 


15.8±0.8 5.3±0.3 




47.5±2.6 


21.2±1.2 7.9±0.7 




0.084 


3.0 


3.0 


0.0 


2005-04-04 


40.5±2.0 


15.2±0.8 5.7±0.3 




46.1±2.5 


21.7±1.2 9.1±0.7 




0.067 


3.0 


3.0 


0.0 


2005-08-18 


42.8±2.1 


14.6±0.7 5.6±0.3 




50.1±2.7 


21.1±1.2 9.2±0.7 




0.080 


3.0 


3.0 
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Observations at 15 GHz 






















1999-11-20 


6.7 ±0.8 


46.2±4.6 




8.1 ±1.3 


50.0±5.0 




0.39 
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2000-03-06 


7.6 ±0.9 


53.6±5.4 




9.5 ±1.5 


55.4±5.6 




0.46 


3.0 


3.0 


0.0 


2000-11-30 


16.2±1.6 


53.7±5.4 




17.6±2.0 


57.8±5.8 




0.27 


3.0 


3.0 
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Observations at 22 GHz 






















2004-11-03 


(28.6) 


(81.9) 




(28.6) 


(80.6) 
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0.0 


2005-04-04 
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0.84 


3.0 


3.0 


0.0 


2005-08-18 


(132.8) 


(62.5) 




(147.1) 


(67.1) 




2.09 


3.0 


3.0 


0.0 
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2000-03-06 


26.41±0.21 


126.2±0.5 


26.17+0.21 


86.9+0.5 


2000-11-30 


26.53±0.21 


126.6±0.5 


25.91+0.21 


87.4+0.5 


2002-09-22 


27.20±0.21 


127.4+0.5 


25.77+0.21 


88.5+0.5 


2004-06-14 


27.77±0.21 


126.9±0.5 


25.89+0.21 


88.7+0.5 


2004-11-03 


27.79±0.21 


126.8+0.5 


26.13+0.21 


88.6+0.5 


2005-04-04 


27.69+0.21 


126.8+0.5 


25.79+0.21 


87.9+0.5 


2005-08-18 


27.85±0.21 


126.5±0.5 


25.79+0.21 


88.5+0.5 


Observations at 15 GHz, 










1999-11-20 


26.82+0.21 


128.0±0.5 






2000-03-06 


26.76±0.21 


127.5±0.5 






2000-11-30 


26.77±0.21 


127.2+0.5 






Observations at 22 GHz 










2004-11-03 


28.24+0.14 


127.0±0.3 






2005-04-04 


28.39+0.14 


127.2+0.3 






2005-08-18 


28.50±0.14 


126.8±0.3 







Table 3. Distances (r) and position angles (PA) of components A and E with respect to B. 




Figure 3. Example of a 5 GHz phase-referenced and self-calibrated im- 
age of NGC 3079, observed on 18 August 2005. At this frequency, all 
components can be seen. Components A, B, and E are still strong, the 
jet-like feature C between A and B is clearly detected, and F is re- 
solved and faint, but detected. Axes are relative to RA 10h01m57.7906s 
Dec 55°40'47.1788" (J2000), and self-calibration produced no measur- 
able shift in the coordinates from phase referencing. Contours are drawn at 
0.24 mjy beam- 1 X 2 N (N = 0, 1, 2, ...). The data have been tapered 
to a resolution of 3 X 3 mas 2 . 



uary 1995 and March 2001. This means that component B has 
been stationary with respect to the maser emission for a decade. If 
the masers trace a molecula r disk around the ce ntral massive black 
hole I Kond ratko etaDl2005l : | Trotter et and their velocity 

dispersion is 400 kms -1 , or 0.0013 c, then masers moving parallel 
to the plane of the sky could have moved by about 0.013 l.y., or 
0.004 pc in 10 yr. This corresponds to 0.06 mas, which would not 
have been noticed. Hence B can be regarded as stationary with re- 
spect to the central black hole powering the AGN, and measuring 
proper motions relative to B is therefore equivalent to measuring 
them relative to the black hole. 



Figure 4. Example of a 15 GHz image of NGC 3079, observed on 30 
November 2000. Components A and B are clearly detected, but there are 
only hints of components C and E. Absolute positions were not measured 
(phase referencing was not used), contours are at 0.81 mjy beam -1 X 
2 N (N = 0,1,2,...). The data have been tapered to a resolution of 
3x3 mas 2 . 

3.3 Separation of A and B 

VLBI observations of NGC 3079 have shown early on that the 
separation of A and B is increasing jMiddelberg et al.l I2005L 
iKondratko et alJ[2005h . We have plotted the previously published 
measurements of the separation of A and B in Fig. [7] together with 
our new measurements. 

The pre-1998 data show a steady expansion of the source, with 
the no table exception of the 5 GHz measurement by iTrotter et al.l 
(1998). The data taken between 1999 and 2003 already indicate 
a slightly slower expansion rate, and the 2004-2005 5 GHz data 
clearly show that the expansion has come to a halt. We note that 
the separation of A and B is consistently larger at higher frequen- 
cies, and we therefore restrict our analysis of the proper motions to 
5 GHz whenever possible. 

The separations measured at higher frequencies do not fol- 
low the same trend as the 5 GHz measurements. The separations at 
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Figure S. Example of a 22 GHz continuum image of NGC 3079, ob- 
served on 18 Aug 2005. Only components A and B have been detected, 
and then with relatively low SNR. Phase referencing to the brightest 
maser did not yield absolute positions at this frequency, contours are at 
6.3 mjy beam -1 X 2^ (N = 0, 1, 2, ...). The data have been tapered 
to a resolution of 3 X 3 mas 2 , which caused a large increase of the image 
noise. The cross indicates the position of the brightest maser feature. 
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Figure 6. A spectrum of the water maser emission in NGC 3079, derived 
from averaging the cross-power spectra of all baselines. No lines were de- 
tected outside the velocity range shown. 
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Figure 7. Upper panel: The evolution of the separation of A and B, us- 
ing data from the literatrure and our measurements. L iterature data are 
denoted with open symbols and have been taken fro m llrwin & Seaquisj 
<198Sl) (epoch 1986.455, 5 G Hz), iTrotter et alj h998l) ( epochs 1992.706, 
5 GHz, and 19 92.745, 8 GHzt.lSawada-Satoh et alJfeOOOl) (epoch 1996.802, 
8 GHz), and iKondratko et al l <2005l) (epochs 1996.919. 5 GHz, and 
2001 .225, 22 GHz) Lower panel: The 5 GHz flux densities of components 
A, B, and E. Only 5 GHz literature data, denoted with open symbols, were 
incorporated in this diagram and were taken from the same publications as 
for the panel above. Measurements of A, B, and E have been connected 
with solid, dashed, and dotted lines to help guide the eye. 



both 15 GHz and 22 GHz are larger than that at 5 GHz; the 15 GHz 
data from 1999-2000 imply a larger separation and do not seem to 
confirm the expansion seen at 5 GHz; the 22 GHz data even seem 
to indicate that the expansion is continuing. However, this appar- 
ent disagreement on the evolution of the separation of A and B 
at low (5 GHz) and high (15 GHz and 22 GHz) frequencies can 
be explained as a consequence of optical depth effects causing a 
frequency-dependent core shift in either or both components. We 
postulate that at 5 GHz we see only an outer layer of the compo- 
nents) where the general kinematic properties of the whole ac- 
tual emission region could be defined. At 15 GHz and 22 GHz, the 
optical depth of the ambient medium allows us to see predomi- 
nantly into a denser emission region with slightly different kine- 
matic properties. 

We have fitted linear functions to all separations observed in 
1985-1996, to our 5 GHz observations from 1999-2002, and to our 
5 GHz observations from 2004-2005. The results are shown in Ta- 



Period 


mas yr 


Av a ng 

mas yr - 1 


P 


A/3 


1985-1996 


0.52 


0.10 


0.12 


0.02 


1999-2002 


0.32 


0.03 


0.08 


0.01 


2004-2005 


0.03 


0.09 


0.01 


0.02 



Table 4. Results of the linear fits to the separation of A with respect to B as 
a function of time, v ang and Av an g are the apparent angular speed and its 
associated error, and /3 and A/3 are the corresponding values of /3 = v/c, 
assuming that the motion is in the plane of the sky at a distance of 15 Mpc. 



ble [4] and are drawn as solid lines in Figure [7] We illustrate the 
motion on the sky of A and E with respect to B in Figure(8] 
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Figure 8. The motions of components A (left panel) and E (right panel) 
with respect to component B. Pluses represent measurements at 5 GHz, 
crosses those at 15 GHz, and squares those at 22 GHz. Arrows indicate the 
direction of motion of a component. The large crosses in the upper left 
corners indicate the errors of the measurements. 



3.4 Motion of component E 

Contrary to the separation of A and B, which appeared to increase 
steadily over more than a decade, the motion of E is more erratic. 
Its separation from component B is plotted in the right panel of 
Figure[8]and an arrow indicates the direction of motion. During our 
observations in 1999 and 2000, E appeared to be moving inwards, 
approximately in direction of the dynamical centre. This relatively 
clear trend stopped between November 2000 and our follow-up ob- 
servations in September 2002. E has since described a figure re- 
sembling a circle, and in our latest measurements in August 2005 
reached a point where it had been in November 2000. 

3.5 Flux densities of A, B, and E 

The evolution of the 5 GHz integrated flux densities of A, B, and 
E is shown in Fig. [7] The flux density of A exhibi ts two sudden in- 
crease s, the first between the 1996 observations of lKondratko et all 
(2005) and our observations in November 1999, and the other be- 
tween our September 2002 and June 2004 observations. The mean 
of the pre-1999 flux densities of A is (7.7 ± 4.1) mjy, of our mea- 
surements between 1999 and 2002 is (26.0 ± 1.6) mjy, and of our 
post-2004 measurements is (48.5 ± 1.9) mjy. More importantly, 
these increases in flux density are coincident with the decelerations 
of component A as described in Section [3~3l The flux density of B 
shows no conclusive evolution over the years, whereas that of com- 
ponent E appears to be increasing slowly. The shape and internal 
structure of E evolved significantly over the years, and part of its 
erratic flux density evolution may be ascribed to this observation. 

We note that in our last 5 GHz epoch, the brightness tempera- 
tures of A, B, and E were 39 x 10 r K, 17 x 10 7 K, and 7 x 10 7 K, 
indicating that the emission is synchrotron emission from AGN or 
supernovae. Brightness temperat ures of HII regions ge nerally reach 
only of the order of 10 4 K (e.g., ICondon et al.| [l991). and are ex- 
tended with sizes of » 0.1" (>> 7.2 pc in NGC 3079). 



3.6 Spectral indices 

We give in Table[5]all spectral indices calculated from the results in 
Table 2, ignoring the 22 GHz flux densities and using the conven- 
tion that S v oc v a . Unfortunately, we do not have data at similar 
frequencies before and after the deceleration events seen in Fig- 
ure |7J so we cannot analyse the spectral evolution of components 
during these events. Although our observations were designed to 
avoid resolution effects, we cannot rule out that the spectral index 



of F is affected by its extent and by some of the 2.3 GHz flux den- 
sity being missing due to the lack of short baselines. 



4 DISCUSSION 

We discuss the observational results as follows. In Section |4~T1 we 
briefly comment on the spectral indices observed in NGC 3079; in 
Section l4~2l we show that components A and B are radio compo- 
nents in a jet which is undergoing interaction with the surround- 
ing material; Section |4~3l demonstrates that the behaviour of A is 
unlikely to be due to jet bending and Doppler boosting; and Sec- 
tion l4.4l shows that components E and F are by-products of inter- 
actions of a jet with the surrounding material. 



4.1 Spectral indices 

During the epochs of 1999 and 2000, the spectral indices ot\f A 
of components B and E remained constant, and component A in- 
creased from -0.90 to -0.27 within the 8.5 months prior to the ob- 
servations in November 2000. 

The spectral indices measured simultaneously in September 
2002 are more interesting. First, A has a strongly inverted spectrum 
with a|;o = 3.37, as has E with 0^3 = 4.49. Second, F has a very 
steep spectrum between 1.7 GHz and 2.3 GHz, and there are several 
similar limits which also deserve explanation. 

The inverted spectra of A and E in the September 2002 epoch 
can be modelled in terms of a free-free absorber in front of the 
components (similarly inverted spectra and limits can be explained 
using the same reasoning, but we lack the data for modelling). In 
the presence of a free-free absorber, the observed flux density, S v , 
of a component can be related to the intrinsic flux density at 1 GHz, 
5*o, of the component by 



S " = S ° \GHz) Xexp (~ T f"' 



(1) 



where ao is the intrinsic spectral index and Tf / is the optical 
depth of the absorber with 



T ff 



■■ 0.0824 



T 



( v 



N+ N- ds 



(2) 



v K J VGHz/ J cm _J cnr pc 

tosterbrockll989l) . The temperature, T, is commonly assumed 
to be 10 4 K, N + and N— are the densities of positive and nega- 
tive charges in the absorber (assumed to be equal), and ds is the 
path length along the line of sight. The integral is called the emis- 
sion measure and is referred to as k hereafter. EquationQ]has three 
unknowns, the intrinsic flux density, So, the intrinsic spectral in- 
dex, ao, and the emission measure, k. It can be solved when the 
flux density has been measured at a minimum of three frequencies. 
When upper and lower limits are ignored, we can solve Equation Q] 
for A, E, and F, taking into account the measurements at 1.7 GHz, 

2.3 GHz, 5.0 GHz, and 15.4 GHz. In the case of A, we use the si- 
multaneous 2.3 GHz and 5.0 GHz measurements from 2002 and 
add the 15.4 GHz measurements from November 2000. Hence vari- 
ability may affect the results. The models are plotted in Figure [9] 

The spectrum of component F is concave and so is not de- 
scribed by a free-free absorbed power law, and we have instead 
fitted a simple power-law to the three spectral points. For com- 
ponent A we find ao = —1.01, So = 305 mjy, and k = 

7.4 x 10 7 pc cm -6 , leading to an optical depth at 5 GHz of 0.83. 
Thus A appears to be a rather normal source of synchrotron emis- 
sion behind a moderately thick free-free absorber. In the case of 
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Epoch SI A B E F 



1999-11-20 


«f 5 °4 


-1.07±0.2 


1.02±0.15 


< 


-1.12±0.12 




2000-03-06 




-0.90±0.2 


0.97±0.15 


< 


-0.75±0.15 




2000-11-30 




-0.27±0.2 


1.00±0.15 


< 


-1.84±0.08 




2002-09-22 


<*h:l 


> 7.80±5.3 






4.49±0.54 


-2.49±0.6 






3.37±0.57 


> 5.22±0.53 




-1.26±0.19 


-0.06±0.34 




»i:Z 


> 4.61±0.06 


> 4.47±0.05 




0.35±0.18 


-0.74±0.24 



Table 5. Spectral indices of NGC3079 components. Limits have been calculated using three times the image rms. 



E, however, the intrinsic spectral index is ao = —4.40, the intrin- 
sic flux density is So = 17 000mJy, and the emission measure 
is 5.3 x 10 7 pc cm -6 , leading to r// at 5 GHz of 0.59. Again the 
properties of the absorber are reasonable, but cto is difficult to ex- 
plain. 

One possible interpretation is that the population of syn- 
chrotron electrons has been cut off from the central engine and 
has aged. The population will lose the high-energy electrons first, 
leading to a sharp drop in the observed spectral energy distri- 
bution towards higher frequencies. The steep spectral index of 
a2,l = —2.49 of F could also arise from such a situation, for 
which we argue in Section |4~4l 

Another possibility is that the electron energy distribution is 
not a power-law, but is a relativistic thermal distribution. Such 
models have been consid ered by, e.g. Jjones & Hardeel dl979h and 
iBeckert & Duschll Jl997h . and they usually reveal an exponential 
cutoff at high frequencies. Both explanations are unsatisfactory be- 
cause they require the electron distribution cutoff to lie very close 
to the frequency at which r// = 1, but we can offer no better solu- 
tion. 

The result is robust against amplitude errors. Only when the 
measurements at 1.7 GHz and 5.0 GHz are increased by as much as 
50 %, and the 2.3 GHz amplitude is decreased by 50 % does Equa- 
tionQ]yield "common" values for a<j of component E. Amplitude 
errors of this magnitude, and in the required sense, are unlikely. 

Another possibility is resolution effects arising from the lack 
of short interferometer spacings and gaps in the (it, w)-coverage, 
both of which could lead to an underestimate of a partially resolved 
emission component. In particular in the case of F we cannot ex- 
clude that the flux is underestimated at 2.3 GHz or 5 GHz. In this 
case the spectral indices a^l and alio given in Table [5] should be 
regarded as lower limits. For example, fitting a simple power-law 
spectrum for F to the three flux densities in Table 2 yields a spec- 
tral index of a — —0.78 ± 0.29, which is close to the canonical 
value of a — —0.7. 

4.2 A jet interacting with the external medium at the location 
of components A and B 

That components A and B are jet features is strongly supported by 
the elongated structure b etween them (comp onent C, Fig. [3] and 
component D, Fig. 3 in iTrotter et al j 1998), by their high speed 
of separation before 2003, and because the line joining compo- 
nents A and B crosses the dynamica l centre of the accretion disk 
in NGC 3079 dKondratko et"ai]|2005l) . This position will be called 
the position of the central engine here after. The spectral proper- 
ties o f A and B lead lTrotter et alj \ 19981) (see also lKondratko et all 
l2005h to suggest that they are lower redshift and lower luminosity 
versions of compact symmetric objects, whose observational prop- 
erties are known to be due to the interaction between jets and their 
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Figure 9. Spectra of components A, E, F, and B, measured in September 
2002 and using 15 GHz data from November 2000. In the cases of com- 
ponents A and E the lines represent the model described in Section |4. 1| 
whereas simple power-laws have been used for B and F. 

high density external medium. These authors reported components 
A and B to have similar spectra, which correspond to synchrotron 
self-absorbed or free-free absorbed regions with t urn-over frequen- 
cies in the range of 8 GHz to 10 GHz. Moreover, Middelber g et al.l 
(2005) reported the spectra of the radio continuum components 
within the innermost 1 pc to 2 pc from the position of the central 
engine to be consistent with an scenario in which the nucleus of 
NGC 3079 is embedded in a dense medium. All these lines of evi- 
dence point to a jet-external medium interaction scenario. 

Our new results show that component B has remained essen- 
tially stationary both in position and in flux density since 2000. 
In contrast, component A underwent two significant flux density 
increases with two corresponding decelerations, the last one com- 
pletely halting the advance of this component. 

Thus, these new results not only agree with the previously 
proposed jet-feature nature of components A and B interacting 
with the dense external medium, they also add further and stronger 
support to this scenario. The "brighter-when-slower" behaviour of 
component A is expected when a propagating jet encounters a pro- 
gressively denser external medium region (e .g. a non-homoge neous 
dense cloud) . The numerical simulations of lWang et al.l j200Ct) and 
ISaxton et al] J2005l) are clear about what happens in that case. If 
the jet is "weak" (i.e. much less denser and/or much less ener- 
getic) with respect to the external dense cloucJE the impact slows 

see IWangetalJbOQd) for the definition of the parameters controlling the 
jet weakness with regard to its external cloud, and for examples of weak 
and powerful jets 
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down the advance of the "hot-spot" of the jet and raises its pres- 
sure. This pressure in crease enhances the synch rotron emission. In 
the approximation of Saxt on et alJ d200lll2002h . an increase of the 
synchrotron emissivit)0by a factor of 10, as experienced by com- 
ponent A between 1997 and 2005 (see Fig. [7J, would occur if the 
pressure increased by a factor of ~3.5. 



4.3 Unlikely jet Doppler boosting scenario for component A 

Note that there is still a scenario that has not been proposed before 
and that might, in principle, explain the "brighter-when-slower" 
behaviour of component A. This scenario also involves a jet but 
the reason for component A's deceleration and brightening is the 
change of propagation direction of component A towards the line 
of sight, together with the Doppler boosting of the synchrotron ra- 
diation. 

If the emission of component A is beamed relativistically then 
the observed flux density, S v , and the intrinsic flux density, S' v i , are 
related by 



(3) 



where S is the Doppler factojfl (e.g. iHughes & Milledll99ll) . 
An increase in the observed flux density of component A by a factor 
of two (as was observed between the end of 2002 and the beginning 
of 2004) would imply an increase of S by a factor of 2 1 / 3 , which 
would also be the minimum required value of the Doppler factor 
after the beginning of 2004. As, by that time the condition 9 ~ 0° 
is required (because the motion of A in the plane of the sky has 
stopped), the intrinsic bulk plasma speed is constrained to be f3 > 
0.23 c. An angle of the jet to the line of sight of ~ 90° (i.e. the jets 
are in the plane of the sky) is usually assumed if the accretion disk 
of the system is almost edge-on, as in the case of NGC 3079. The 
large acceleration required both to bring the bulk plasma from j3 ~ 
0.13 c to f3 « 0.23 c and to bend the jet beam by an amount of 45° 
or larger, make the jet-Doppler boosting scenario rather unlikely. A 
lower angle (> 45°) is assumed here to account for the fact that the 
jet is cl early not aligned with the angular momentum of the disk 
(see e.g. lKondratko et alj2005h but at ~ 45° from such direction. 

Hence, we conclude that the Doppler boosting scenario re- 
quires the jet in NGC 3079 at the location of A to bend by ~ 45° 
in a time range of only one year. In such conditions, a strong im- 
pact with the external medium unavoidably produces such extreme 
bending. In this case, the Doppler boosting scenario brings us again 
to the interaction of the jet with the external medium. 



4.4 Components E and F as by-products of jet-external 
medium interaction 

The fact that components E and F were found at positions 
not aligned wit h the j et between components A and B made 
iKondratko et alJ ( 120051) suggest that E and F were intrinsically 
different from A and B. In addition, the spectral properties of E 
and F (with much steeper spectra t han those correspondin g to op- 
tically thin synchrotron spectra) led lKondratko et aU d2005) to pro- 



2 j v oc p( 3 a )/ 2 , where p is the pressure and a is the spectral index of 
the synchrotron emission. We assume a ss —0.7 here. 

3 5 = [T(l - /3COS9)]- 1 where F = (1 - /J 2 )" 1 / 2 is the Lorentz factor, 
P is the speed in units of c and 8 is the angle between the direction of the 
flow and the line of sight. 



pose that E and F are consistent with ageing synchrotron sources; 
i.e. with no injection of fres h non-thermal electrons. 

IKondratko etalj ( [20051 showed that the possibility that E and 
F are radio supernovae is highly unlikely. Instead, they proposed 
two different scenarios to explain the properties of E and F: a) ei- 
ther they are remains of emission from a wobbling jet which was 
ejected earlier in orientations different from the current jet direc- 
tion or b) they are shocks in a wide-angle, pc-scale outflow inter- 
acting with its external medium. Within the latter picture, the out- 
flow would be a pc-scaled version of the kpc-sc ale bipolar super- 
bubbles observed in radio c ontinuum emission touric & Se aquist 
1 19881 ; iBaan & Irwml 1 19951; lirwin & Saikial 120031) , in TN III+Hq 
emission dVeilleux et al.ll 19941; ICecil et al.ll200ll) and in soft X-ray 
continuum emission ( Cecil et al.| [2002). 

We consider the jet wobbling scenario (a) very unlikely since 
no evidence has been found to support a systematic change of the 
position angle of the current jet. In this scenario, the jet axis defined 
as the line which connects components A and B, should systemati- 
cally change with time, which has not been observed between 1994 
and 2005 (cf. Table [3}. Furthermore, component C nicely line s up 
with the A-B axis in all images since lirwin & Seaquistl Jl988h ob- 
served NGC 3079 for the first time with VLBI. Hence, during this 
time range, the jet wobbling scenario can be safely ruled out. 

Although scenario (b) is not as unlikely as the jet wobbling 
scenario, we still find it difficult to imagine an outflow ejected from 
the central engine in NGC 3079 which interacts with the dense ex- 
ternal medium and produces a shock moving inwards, in the oppo- 
site sense to the assumed outflow, as detected for component E in 
1999-2000 (Figure[8). 

Here, we propose an alternative scenario to account for the 
presence of synchrotron emission regions from non-thermal aging 
electrons close to the pc-scale jets in NGC 3079. This new scenario 
comes up naturally from the jet/dense-external-medium interaction 
picture previously proposed to explain the behaviour of component 
A. Our idea is based on the two-dimensional numerical simulatio ns 
of jets in inhomogeneous media reported by Saxton et 

]d]|2005fl 

Their grid of simulations, on which a tenuous medium is also 
populated by a relatively small filling factor of dense clouds, show 
that when the filling factor of the clouds is large, the jet has more 
possibilities to be disrupted and to form channels through which the 
shocked radio emitting plasma can flow between the clouds. The 
formation of a hierarchical structure of nested lobes (not always 
producing strong radio emission) and an outer enclosing bubble 
seems unavoidable in all the simulated conditions. Having in mind 
this scenario, it is difficult not to find a close, at least qualitative, 
relation between our component A with the region of the simulated 
jet in which the first dense cloud is encountered and the jet is dis- 
rupted; and between components E and F with the simulated lobes 
or knots of already shocked material. In the simulations, this ma- 
terial slowly flows from the first shocked regio n along the opened 
channels between the clouds (see e.g. Fig. 1 1 of lSaxton et al.l2005l) 
and not necessarily outwards, away from the central engine, be- 
cause their motions depend on the cloud-obstacles that they en- 
counter in propagation. This could explain why component E does 
not move outwards but in wards, towards the central engine (Fig- 
ure [8). The simulations of lSaxton et alJ d2005l) show that the lobes 
of shocked material can be found at large distances (relative to the 



4 The scaling of the results from Saxton et al. ( 2005) was set to apply their 
quantitative results to kpc-scale jets. However, for us who use their study 
on a qualitative way only, their results are equally applicable. 
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jet size) from the first shocked region and in directions that are far 
from where the original straight jet would have propagated. 

A consequence of this is that, as the clouds move, plasma 
channels that initially emit synchrotron emission can eventually be 
closed, cutting off synchrotron-emitting regions from a fresh sup- 
ply of non-thermal particles. As a consequence, the electron pop- 
ulation of synchrotron-emitting regions is expected to age. This 
process manifests itself in an overall decrease of flux density and 
the steepening of radio spectra, as high-energy particles loose their 
energy more quickly and radio emission at higher frequencies de- 
creases first. In NGC 3079, the 5 GHz flux density of component E 
has approximately doubled between 1999 and 2005, which seems 
to disagree with the scenario of an aging electron population. How- 
ever, we note that E underwent significant internal evolution in 
1999 and 2000. This evolution maybe caused by compression of 
the component, which caused the observed increase in flux density. 



5 SUMMARY AND CONCLUSIONS 

Based on the similarities of the spectra of GPS sources with those 
of components A and B, previous studies proposed these com- 
ponents to be regions in which a jet, launched from the central 
super-massive engine i n NGC 3079, interacts w ith the dense exter- 
nal medium. Moreover, Middelb erg et alj i2005h reported the spec- 
tra of the radio continuum components within the innermost 1 pc 
to 2 pc from the accretion disk to be consistent with an scenario in 
which the nucleus of NGC 3079 is embedded in a dense medium. 
We have shown that the structure and the time evolution of the ra- 
dio continuum emission in the innermost pc around the accretion 
disk in NGC 3079 behaves as expected for a jet interacting with 
a clumpy external medium. Hence, our new results not only agree 
with both the jet scenario and the jet-interaction scenario previously 
proposed, but also strongly support both of them. 

Within our proposed scenario, depending on the properties of 
the cloudy medium in which they are embedded, A and B are ex- 
pected to continue separating from the dynamic centre of the accre- 
tion disk once they are able to drill through the cloud that stopped 
them, as was proposed for the expanding phase of the jet in III Zw 2 
jBrunthaler et al.l2005l) . Note however, that the re-activation of the 
expansion of components A and B might start to be significant after 
large time periods (depending on the relative properties of the jets 
and the clumpy medium on which they are embedded). Thus, these 
expansions may not be observable for years or tens of years. The 
spectra of components E and F should show evidence of a slowly 
ageing non-thermal electron population until they have faded com- 
pletely. Although not required, they can be expected to move with 
speeds smaller than those previously measured for A and in di- 
rections which do not necessarily have to point outwards from the 
central engine in NGC 3 079. 

The simulations of ISaxton et ail J2005h show that, if the fill- 
ing factor of dense clouds in the surrounding medium of the jet is 
high, the momentum of the jet spreads at the first strong interac- 
tion and forms a slowly expanding bubble which grows larger than 
the jet itself, which indicates that the jet is impeded from propa- 
gating freely. If one assumes a large filling factor of dense clouds 
surrounding the central region of NGC 3079, the interaction sce- 
nario offers, as a natural consequence of the slower propagation 
of the jet, the explanation for the origin of the kpc-scale bipolar 
super-bubbles observed in radio continuum emission, [N II]+Ha 
emission and soft X-ray continuum emission in NGC 3079. 

There is at present little doubt about the presence of strong 



interactions of th e radio plasma with th e inter-stellar medium in 
Seyfert galaxies dOosterloo etal] l200d) . As for NGC 3079, ev- 
idence for interaction of the jets with their external medium 
in subgalactic-scales seems t o be common in Seyfert galaxies 
(e.g. NGC 1068 and N GC 4 1 5 l. lWilson & Ulvestad 19821 ; III Zw 2, 
iBrunthaler et alj [2005). Hence, it appears reasonable to apply the 
interaction scenario to the whole class of Seyfert galaxies. This 
would give a satisfactory explanation for why in a large fraction 
of Sey fert galaxies large-scale bipolar super -bubbles are observed 
(e.g. IColbert et alJI 19961 IColbert et al.ll 19961) instead of well colli- 
mated, kpc-scale jets as in radio-loud AGN. 
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